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1.  INTRODUCTION 

Breaking  waves  are  a  ubiquitous  phenomenon 
of  the  world’s  oceans.  They  disrupt  the  aqueous 
boundary  layer  causing  surface  renewal,  thereby 
enhancing  the  diffusion  of  gases  and  heat  across 
the  air-sea  interface.  Breaking  waves  are  also 
responsible  for  the  dissipation  of  wave  energy  and 
thus  directly  affect  the  evolution  of  the  wind-wave 
spectrum.  With  advances  in  technology,  new  direct 
observations  of  the  two-dimensional  spatial  surface 
wave  topography  have  been  made.  These  data 
allow  for  the  opportunity  to  go  beyond  linear 
analysis  and  study  the  nonlinearity  of  the  surface 
wave  field,  in  particular  the  statistics  of  steep  and 
breaking  waves. 

While  most  previous  observational  studies  have 
used  whitecaps  or  bubbles  to  detect  breaking  wave 
events,  new  mathematical  methods  of  data  analysis 
have  been  exploited  in  trying  to  detect  breaking 
waves  from  wave  height  records.  In  particular,  Liu 
(1994)  used  the  Morlet  wavelet  transform  to 
analyze  ocean  wave  data  taken  with  a  wave  wire 
mounted  on  a  buoy.  Although  Uu’s  (1994) 
approach  is  limited  to  a  narrow  band  wave  system 
and  cannot  be  applied  to  open  ocean  surface  wave 
data,  the  wavelet  transform  approach,  in  principle, 
should  be  suitable  for  detecting  wave  breaking. 

This  is  because  previous  theoretical  studies  (Dold 
and  Peregrine  1986;  Banner  and  Tian  1998) 
suggest  that  breaking  wave  events  are  associated 
with  wave  groups  with  strong  nonlinearity  rather 
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than  with  a  single  steep  wave,  and  wave  groups 
can  be  detected  by  the  Morlet  wavelet  transform.  In 
this  study  we  propose  a  new  approach  to  estimate 
the  statistics  of  steep  wave  events  (wave  groups  of 
large  amplitude)  by  applying  the  Morlet  wavelet 
transform  to  broad-band  open  ocean  wave  fields. 
We  make  use  of  the  spatial  wave  topographic  data 
obtained  during  the  Southern  Ocean  Wave 
Experiment  (SOWEX)  in  June  1992  and  the 
experiment  conducted  off  Duck,  North  Carolina  in 
September  1997  (hereafter  termed  the  DNC 
experiment).  These  data  were  obtained  with  a 
Scanning  Radar  Altimeter  and  an  airborne  scanning 
lidar  system  respectively  which  produce  two- 
dimensional  topographic  maps  of  the  sea  surface 
wave  height  (see  Banner  et  al.  1999  and  Hwang  et 
al.  2000).  The  results  are  then  used  to  examine 
how  the  statistics  of  nonlinear  wave  groups  may 
correlate  with  the  true  breaking  wave  statistic.  In 
this  study  we  present  a  one-dimensional  data 
analysis  based  on  the  assumption  that  all  steep 
waves  propagate  in  the  mean  wind  direction. 

2.  DATA  ANALYSIS 

2.1  Definition  of  the  steep  wave  statistic 

Phillips  (1985)  proposed  a  distribution  function 
A (c)  such  that  A(c)dc  represents  the  average 
total  length  per  unit  surface  area  of  breaking  wave 
fronts  that  have  intrinsic  velocities  in  the  range  c  to 
c  +  dc  .  In  the  same  manner,  we  define  the  steep 
wave  statistic  Ar  (k)  as  the  total  length  of  steep 
wave  fronts  with  wave  slope  exceeding  a  set 
threshold  T  per  unit  surface  area  per  unit 
wavenumber.  The  steep  wave  statistic,  A T(k)  is  a 
dimensionless  quantity  which,  when  integrated  over 
all  wavenumbers  yields  F ,  the  total  length  of  steep 


20060130  276 


wave  fronts  above  wave  slope  threshold  T  per  unit 
surface  area.  In  general,  the  slope  of  an  individual 
wave  cannot  be  determined  uniquely  for  random 
seas  with  a  broad-banded  spectrum.  Thus,  the 
wave  slope  is  defined  using  the  wavelet  transform 
such  that  the  estimated  slope  of  an  individual  wave 
is  in  fact  the  average  wave  slope  of  a  small  group 
of  waves  that  are  detected  by  the  wavelet 
transform. 

The  wavelet  transform  of  signal  f(u )  in  this 
study  is  defined  as 


Wf(a,s)  =  ' W  a >0(1) 
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is  the  Morlet  wavelet,  a  is  the  scale,  and  s  is  the 


underwent  the  wavelet  transform.  Wavelet  analysis 
was  performed  over  a  set  of  discrete  scales  (with  a 
constant  A  a )  that  extend  from  a  scale  larger  than 
the  dominant  wavelength  down  to  a  scale  before 
the  Fourier  spectral  noise  floor  in  the  data  appears. 

Each  column  vector  of  the  wave  height  image 
admits  a  two-dimensional  function,  Wf(a,s ) 
dubbed  the  wavelet  transform.  There  is  a  one  to 
one  relationship  between  the  wavenumber  and  the 
wavelet  scales.  The  wavenumber  associated  with 
the  Morlet  wavelet  at  scale  a  is  taken  to  be  the 
wavenumber  associated  with  the  peak  in  the  power 
spectrum  of  the  Morlet  wavelet.  The  conversion 
from  the  wavelet  transform  threshold  to  the 
equivalent  local  wave  slope  threshold  is 
accomplished  by  first  taking  a  pure  sinusoid 
f(u)  =  B  sin  ku  ,  of  known  wave  amplitude,  B , 
and  wavenumber,  k  and  obtaining  the  wavelet 
transform  of  the  signal,  Wf(a,s).  The  peak  value 
of  the  wavelet  transform  is  then  related  to  the  wave 
slope  Bk  ,  to  determine  the  conversion  constant 
ry. 


location.  Note  that  the  classical  _ ...  normalization 
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is  not  used.  The  Morlet  wavelet  transform 
possesses  two  aspects.  When  applied  to  a  signal  it 
effectively  searches  the  signal  and  finds  regions 
where  the  data  has  wave-like  characteristics.  In 
addition,  if  the  signal  is  interpreted  as  wave  height, 
the  local  peak  value  of  the  inner  product  of  the 
signal  and  the  Morlet  wavelet  produces  a  measure 
of  the  average  wave  slope  over  the  support  of  the 
wavelet.  Thus,  a  signal  with  a  high  wave  slope 

event  with  its  scale  a0  and  its  crest  position  s0  will 
have  a  wavelet  transform  characterized  by  a  large 
peak  value  of  Wf  at  (a. ) ,  and  the  peak  value  is 
proportional  to  the  average  slope  of  the  signal  in 
the  neighborhood  of  j  . 

2.2  Estimation  of  the  steep  wave  statistic, 

Ar(*)- 

2.2.1  Wavelet  analysis  of  signal 

Initially,  the  SOWEX  and  the  DNC  data  sets 
were  preprocessed  for  the  purpose  of  eliminating 
high  frequency  noise,  and  dropouts.  The  SOWEX 
data  was  divided  into  subsets  according  to 
experimental  day.  This  data  along  with  the  DNC 
data  set  were  rotated  and  bilinearly  interpolated 
such  that  the  y-axis  is  in  the  mean  downwind 
direction  and  such  that  all  images  of  an  individual 
subset  possess  the  same  resolution  in  the  along- 
wind  direction  and  cross-wind  direction.  For  each 
data  set  analyzed,  columns  of  the  data,  which 
constitute  data  in  the  direction  of  the  wind,  were 
stripped  off.  These  data  vectors  subsequently 
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A  contour  plot  of  the  wavelet  transform  along 
with  a  plot  of  an  analyzed  signal  is  shown  from  the 
June  10th  SOWEX  experiment  in  Figure  1. 


Figure  1:  June  10th  SOWEX  experiment  surface 
elevation  signal  (bottom  panel),  and  contour  plots  of 
its  wavelet  transform  (top  and  middle  panels).  Solid 
contour  lines  are  drawn  at  every  0.04.  Dashed 
contour  lines  are  at  wave  slope  threshold  of  0.02 
(top  panel)  and  0.06  (middle  panel).  Wave  crests 
exceeding  wave  slope  threshold  of  0.02  (top  panel) 
and  0.06  (middle  panel)  are  marked  with  asterisks. 


The  wavelet  transform  Wf(a,s )  contains  an 
array  of  high  wave  slope  events  which  appear  as 
local  maxima  throughout  a,s  space.  ln,order 


to  obtain  a  distribution  of  these  events,  a  wave 
slope  threshold  is  applied  over  Wf(a,s ) 
designated  by  the  dashed  lines.  The  local  maxima 
that  satisfy  the  condition  of  being  above  the  set 
threshold  appear  in  the  figure  as  asterisks.  These 
asterisks  are  defined  as  individual  steep  wave 
events  in  the  subsequent  analyses.  It  is  worthy  of 
note  that  more  than  one  local  maximum  (asterisk) 
can  exist  within  a  single  region  enclosed  by  a 
dashed  line  when  the  threshold  is  set  low. 

Each  high  wave  slope  event  in  a  data  vector, 
denoted  by  asterisks  in  Figure  1,  has  associated 
with  it  a  scale  value  and  an  along  wind  position 
value.  The  scale  values  together  form  DT  (a) , 


where  DT  is  the  number  of  events  of  scale  a 

exceeding  a  set  wave  slope  threshold,  T  .  The 
same  method  of  analysis  outlined  above  is  applied 
to  all  of  the  data  vectors  (with  index  n)  that  are  part 
of  the  full  two-dimensional  wave  height  image  to 
obtain  a  set  of  distributions,  D”  (a) .  This 
distribution  is  easily  transformed  into  the  lambda 
function,  A r  (a)  by  summing  the  distributions  over 
n  and  then  normalizing  the  result.  Thus,  the  final 
result  for  Ar  (a)  is 


A  r(a) 


N 

^Dr(a)x  Ax 

n= 1 


AxAa 


(4) 


Here,  the  cross-wind  sampling  distance,  Ax  is 
taken  to  be  the  length  of  the  high  wave  slope  crest, 
A  a  is  the  differential  scale,  N  is  the  total  number 
of  data  vectors,  and  A  is  the  total  area  of  the  image 
analyzed.  The  lambda  function  is  averaged  over  a 
set  of  realizations  and  then  finally  transformed  into 

Ar(*). 


2.2.2  Doppler  correction  for  the  steep  wave 
statistic 

The  Doppler  shift  of  the  waves  due  to  the  aircraft 
motion  is  a  major  source  of  error.  The  error 
associated  with  this  effect  is  estimated  in  terms  of 
the  apparent  and  actual  value  for  Ar  (k) .  This 
Doppler  correction  was  performed  assuming  that  all 
wave  fronts  are  perpendicular  to  the  mean  wind. 

3.  RESULTS  AND  DISCUSSION 

3.1  Detection  of  steep  wave  crests. 

The  wavelet  analysis  algorithm  of  the  surface 
wave  topography  data  vectors  is  successful  at 
pinpointing  high  wave  slope  crests.  Figure  1  shows 


a  space  series  taken  from  June  10th  of  the  SOWEX 
experiment  along  with  a  contour  plot  of  the  wavelet 
transform  values.  Asterisks  identify  high  wave  slope 
events  that  exceed  a  wave  slope  threshold  of  0.02 
(top  panel)  and  0.06  (middle  panel).  The  algorithm 
is  able  to  perform  a  multi-scale  filtration  of  the  data 
to  reveal  steep  wave  crests  at  both  large  and  small 
scales.  The  steep  wave  slope  events  identified  by 
the  asterisks  at  small  and  large  wavenumbers  are 
clearly  correlated  with  large  and  small  steep  waves 
respectively  in  the  surface  height  signal.  Some 
asterisks  do  not  lie  directly  at  the  points  of 
maximum  surface  elevation  in  the  signal.  This  is 
due  to  the  inner  product  associated  with  the 
convolution.  It  attempts  to  make  the  best  fit  of  the 
Morlet  wavelet  to  the  data  and  occasionally  falls 
short  of  locating  wave  crests  in  agreement  with  the 
actual  crest  locations. 

3.2  Characteristics  of  the  steep  wave  statistic 

The  estimates  of  the  steep  wave  statistic  with 
the  Doppler  correction  ,  AcT  (k)  for  three  SOWEX 
cases  and  the  DNC  experiment  are  shown  in  Figure 
2.  The  case  June  12th  is  not  shown  here  since  it  is 
similar  to  June  10th. 
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Figure  2:  Corrected  steep  wave  statistic  AcT  (k) 
vs.  k  for  different  wave  slope  threshold  values. 

Numbers  in  the  figure  indicate  slope  threshold  T . 
(a)  DNC  experiment.  Mean  wind  speed  9.5  m/s.  (b) 
June  10th  SOWEX  experiment.  Mean  wind  speed 
20  m/s.  (c)  June  13th  SOWEX  experiment.  Mean 
wind  speed  9  m/s.  (d)  June  14th  SOWEX 
experiment.  Mean  wind  speed  =  6  m/s. 

The  curves  in  Figure  2  in  general  display  low 
values  of  AcT  ( k )  below  the  dominant  wavenumber 
and  increased  amounts  of  steep  wave  crests  at 
high  wavenumbers.  These  results  bear  qualitative 
resemblance  to  the  breaking  wave  results  of  Ding 
and  Farmer  (1994).  They  found  that  the  average 
scale  of  breaking  to  be  much  smaller  than  the 
dominant  wind  wave.  The  steep  wave  statistic  is 
consistent  with  this  result  with  the  peak  in  the  steep 
wave  statistic  occurring  at  a  much  higher 
wavenumber  than  the  dominant  wind  wave.  What  is 
clearly  shown  in  these  plots  is  that  all  data  sets 
have  approximately  the  same  average  length  of 
wave  crests  per  unit  area  per  unit  wavenumber, 
regardless  of  the  wind  forcing  or  the  wave  field,  at 
the  lowest  wave  slope  threshold  of  0.01.  Thus, 
statistically  the  total  number  of  wave  crests  is 
similar  for  all  cases  if  you  count  all  wave  crests  (low 
and  high  wave  slope  waves)  in  the  wave  field. 

Differences  in  AcT(k)  become  apparent  with 
increases  in  the  wave  slope  threshold  T.  For  the 
SOWEX  data  set,  Figure  2  clearly  shows  that 
AcT  (/t)  at  various  wavenumbers  decreases  much 


more  rapidly  with  increasing  T  for  the  low  wind 
cases  than  for  the  high  wind  cases.  This  is  to  be 
expected  since,  to  the  first  order,  the  higher  wind 
cases  tend  to  possess  larger  amounts  of  high  wave 
slope  events.  Figures  3  a)  and  3  b)  show  that 

AcT  (k)  at  a  fixed  k  from  both  the  SOWEX  and 
DNC  data  sets  decays  exponentially  with  wave 
slope  threshold  squared,  i.e.,  it  approximately 
follows  the  form  of 

AcT  =  Aae~pT'  (5) 

with  a  fixed  exponent  p  and  multiplicative 
constant  of  Aa .  We  have  made  a  fit  of  the  data  to 
this  function  and  chose  AcT  (k  =  0.04  rad  /  m) 
from  the  data  of  the  SOWEX  experiment  and 
A cT(k  =  0.4  radhri)  from  the  DNC  data,  since 
these  wavenumbers  are  significantly  higher  than 
the  dominant  wavenumber,  but  above  the  point  of 
noise  for  all  the  data  sets. 


a) 


b) 

Figure  3:  The  steep  wave  statistic  AcT  (k)  vs. 

wave  slope  threshold  squared  ( T 2 ).  Data 
designated  by  solid  lines.  Least  square  fit 
designated  by  dots  connected  by  solid  lines,  (a) 
SOWEX  experiment.  k  =  0.04  rad  lm  .  (b)  DNC 
experiment.  k  =  0Arad/m.  Thick  solid  line  is 
Phillips  et  al.'s  (2001)  estimate  of  A .  Thick  dashed 
line  is  Melville  et  al.’s  (2002)  estimate  of  A . 


3.3  Relationship  between  the  steep  wave 

statistic  and  the  breaking  wave  statistic 

If  it  is  conjectured  that  the  steep  wave  statistic  at 
a  particular  high  wave  slope  threshold  is  equivalent 
to  the  true  breaking  wave  statistic,  the  trend  of 
AcT  (£)  from  the  SOWEX  data  sets  implies  that 
there  are  more  breaking  waves  at  higher  wind 
forcing  conditions.  Furthermore,  if  the  breaking 
wave  statistic  A(k)  is  known  independent  of 

it  is  possible  to  estimate  a  breaking  wave 
slope  threshold  by  comparing  the  two.  Phillips  et  al. 
(2001)  calculated  A(e),  with  a  correction  for  the 
advection  of  breaking  wave  crests  by  swell,  over  a 
range  of  wave  speeds  from  2.5  m/s  to  6  m/s  at  the 
mean  wind  speed  of  9.3  m/s.  Melville  et  al.  (2002) 
reported  A(c)  at  three  wind  speeds  of  7.2,  9.8, 
and  13.6  m/s.  Therefore,  the  results  of  Phillips  et  al. 
(2001),  Melville  et  al.  (2002)  (case  of  wind  speed 
9.8  m/s),  and  the  DNC  data  set  all  have  a  similar 
wavenumber  range  and  were  obtained  under 
approximately  the  same  wind  forcing  conditions  of 
Ul0  =9  mis.  Figure  3  b)  shows  a  linear  fit  to  the 
natural  logarithm  of  AcT  (k-k0)  plotted  with  the 
square  of  the  wave  slope  threshold  at  k0  =  0.4 
rad/m  for  the  DNC  data  set.  The  horizontal  dashed 
and  solid  line  represents  the  value  of  the  breaking 
wave  statistic  at  k0  =  0.4  rad/m  obtained  by 
Melville  et  al.  (2002)  and  Phillips  et  al.  (2001).  The 
intersection  between  the  linear  fit  and  the  horizontal 
lines  gives  two  estimates  of  the  wave  slope 
threshold  at  which  the  steep  wave  statistic 
becomes  identical  to  the  breaking  wave  statistic. 

The  corresponding  wave  slope  threshold  is 
between  approximately  0.11  and  0.13,  with  0.12 
being  the  mid-value.  These  values  are  much  lower 
than  the  traditional  wave  slope  threshold  of  0.4  first 
proposed  by  Stokes. 

Dold  and  Peregrine  (1986)  conducted  numerical 
experiments  in  which  two-dimensional  nonlinear 
wave  packets  of  different  initial  wave  slope  were 
allowed  to  propagate.  They  were  able  to  obtain  the 
result  of  the  wave  slope  threshold  above  which 
wave  groups  developed  into  breaking  waves  versus 
the  number  of  waves  in  a  wave  group.  It  is  of 
interest  to  note  that  for  a  wave  packet  containing  5 
waves,  the  corresponding  wave  slope  threshold  is 
approximately  0.11  according  to  Dold  and 
Peregrine  (1986).  This  value  is  dose  to  our 
breaking  wave  slope  threshold  extrapolated  from 
the  data. 

Assuming  that  the  steep  wave  statistic  at  the 
equivalent  slope  of  0. 12  agrees  with  the  breaking 
statistic,  we  may  estimate  the  breaking  statistic  for 
all  the  SOWEX  experimental  cases.  Figure  4 
shows  the  estimates  of  AcT  (k)  for  the  DNC 
experiment  and  SOWEX  experiment  at 


Figure  4:  Estimates  of  Acr  (A:)  multiplied  by 

(lO/£/10)3  at  wave  slope  threshold  r  =  0.12  for 
the  SOWEX  and  DNC  experiment  (thick  lines). 
Phillips  et  al.’s  (2001)  estimate  of  A(£)  are 
designated  by  crosses.  Melville  et  al.’s  (2002) 
estimates  of  A(&)are  shown  at  U]0  =7.2  m/s 

(-),  9.8  m/s  (•  -)  and  13.6  m/s  (*  -) . 

The  SOWEX  results  of  June  14th  are  not  shown 
since  this  estimate  is  contaminated  with  noise.  Also 
plotted  are  Phillips  et  al.’s  (2001)  and  Melville  et 
al.’s  (2002)  estimates  of  the  breaking  wave  statistic. 
In  order  to  remove  the  wind  speed  dependence, 
both  AcT  (£)  and  A(k)  are  scaled  with  the  cube  of 
the  wind  speed,  as  suggested  by  Melville  et  al. 
(2002).  The  results  of  the  DNC  experiment  are 
close  to  that  of  Phillips  et  al.’s  (2001)  and  Melville  et 
al.’s  (2002)  result  at  Ul0  =  9.Sm/s  as  expected. 
The  SOWEX  results  of  June  10th  and  12th  are  much 
higher  because  the  wind  speed  was  much  higher 
on  these  two  days.  All  results  appear  to  collapse 
within  one  order  of  magnitude  or  so.  This  suggests 
that  the  number  of  steep  waves  at  T  =  0. 12  also 
scales  roughly  with  the  cube  of  the  wind  speed. 
Nevertheless,  some  variability  still  remains  which 
most  likely  is  the  result  of  the  wave  fields  being  at 
different  stages  of  development  on  each  day. 

4.  CONCLUSIONS 

The  wavelet  analysis  methodology  presented 
here  is  able  to  detect  steep  wave  events  and  give 
estimates  of  the  amount  of  high  wave  slope  events 
that  cover  a  given  area  of  ocean.  Analysis  of  the 
results  shows  that  high  wave  slope  crests  appear 
over  the  entire  range  of  wavenumbers  resolved, 
with  a  large  amount  being  much  shorter  in 
wavelength  than  the  dominant  wave.  At  low  wave 
slope  thresholds,  the  total  crest  length  is 


approximately  independent  of  wind  forcing  for  all 
wave  fields  considered.  The  steep  wave  statistic 
AcT  (£)  then  decays  exponentially  with  the  square 
of  the  slope  threshold  T.  If  the  steep  wave  statistic 
is  hypothesized  to  evolve  into  the  breaking  wave 
statistic  at  a  specific  wave  slope  threshold, 
comparison  of  Ac7.(£)with  previous  independent 
measurements  of  the  breaking  wave  statistic  gives 
a  wave  slope  threshold  of  about  0.12.  This 
threshold  is  consistent  with  the  results  of  the 
numerical  studies  of  Dold  and  Peregrine  (1986). 
Comparison  of  the  steep  wave  statistic  at  this 
extrapolated  wave  slope  threshold  with 
independent  breaking  wave  measurements  suggest 
that  breaking  and  steep  wave  statistics  are  both 
roughly  proportional  to  the  cube  of  the  wind  speed, 
but  other  factors  besides  the  wind  speed  also  affect 
the  level  of  the  steep  wave  statistic. 
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